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Aiming at the investigation of above-threshold ionization in super-strong laser fields with highly charged ions, 
we develop a Coulomb-corrected strong field approximation (SFA). The influence of the Coulomb potential 
of the atomic core on the ionized electron dynamics in the continuum is taken into account via the eikonal 
approximation, treating the Coulomb potential perturbatively in the phase of the quasi-classical wave function 
of the continuum electron. In this paper the formalism of the Coulomb-corrected SFA for the nonrelativistic 
regime is discussed employing velocity and length gauge. Direct ionization of a hydrogen-like system in a 
strong linearly polarized laser field is considered. The relation of the results in the different gauges to the 
Perelomov-Popov-Terent'ev imaginary-time method is discussed. 
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I. INTRODUCTION 

Due to advances in laser technology strong near-infrared 
laser fields nowadays are available up to intensities of 
10 22 W/cm 2 |1] and much stronger laser fields are envisaged in 
near future [2] stimulating the investigation of the relativistic 
regime of laser-atom interaction in ultra-strong fields. The pio- 
neering experiment in this field was carried out by Moore et al. 
|3|. They have investigated the ionization behavior of atoms 
and ions in a strong laser field at an intensity of 3 x 10 18 W/cm 2 . 
Several further experiments have been devoted to relativistic 
laser-induced ionization PH- fiOl . 

Numerical investigation of the dynamics of highly-charged 
ions in super-strong fields has been carried out in IfTTlETl . The 
standard analytical approaches in the field of nonperturbative 
laser-atom interaction are the strong field approximation (SFA) 
Il22l - l24l and the imaginary time method (ITM) [25 26 1. For 
a theoretical treatment of the relativistic effects, the SFA has 
been generalized into the relativistic regime in l27l |2§1 and 
the ITM in I29H321 , respectively. In the standard SFA the 
influence of the Coulomb field of the atomic core is neglected 
in the electron continuum dynamics and the latter is described 
by the Volkov wave function 11331 . Accordingly, the predictive 
power of the SFA is the best for negative ions where no long 
range forces of the parent system act on the ionized electron. 
For atoms or molecules with long range Coulomb forces the 
performance of the SFA downgrades to a qualitative level l34l . 
This is true especially for highly-charged ions. 

In the non-relativistic regime the ITM has been successfully 
used to treat Coulomb field effects during the ionization in the 
quasi-static regime and the well-known quantitatively correct 
Perelomov-Popov-Terent'ev (PPT) ionization rate has been 
derived Il35ll36l . The PPT theory uses the quasi-classical wave 
function for the description of the tunneling part of the elec- 
tron wave packet through the quasi-static barrier formed by 
the laser and atomic field, with matching of the quasi-classical 
wave function to the exact bound state wave function 0711381 . 
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The standard SFA technique has also been modified to include 
Coulomb field effects of the atomic core. The simplest heuris- 
tic approach is the, so-called, Coulomb- Volkov ansatz in which 
the Volkov wave function in the SFA matrix element is replaced 
by an heuristic Coulomb-Volkov wave function 1 39-49]. In the 
latter the Coulomb field is taken into account via an incorpora- 
tion of the asymptotic phase of the exact Coulomb-continuum 
wave function into the phase of the heuristic Coulomb-Volkov 
wave function llBTjI . Consequently, the coupling between the 
Coulomb and laser field is neglected in the Coulomb-Volkov 
ansatz and the approach fails when the electron appears in the 
continuum after tunneling close to the atomic core BD . 

Following a more rigorous approach, the eikonal approxima- 
tion ll52l has been proposed to apply for strong field problems 
ll53l . In the latter, nonrelativistic free-free transitions in the 
laser and the Coulomb field have been considered employing 
an eikonal wave function for the continuum electron. Here the 
laser field is taken into account exactly, while the Coulomb 
field is via the eikonal approximation. The eikonal approxima- 
tion has been generalized in [54] to include quantum recoil ef- 
fects at photon emission and absorption. A Coulomb-corrected 
SFA for nonrelativistic ionization employing the eikonal wave- 
function has been first proposed in [55 1. Similar approaches 
have been considered in 1 56-61 1. Recently, the nonrelativistic 
Coulomb-corrected SFA based on the eikonal- Volkov wave 
function for the continuum electron has been further elaborated 
in 1621 163 1 and applied for molecular strong field ionization and 
high-order harmonic generation. The Coulomb-corrected SFA 
has also been extended to include rescattering effects If64ll65ll . 
Here the Coulomb field is taken into account exactly in the 
quasi-classical electron continuum trajectories that are later 
plugged into the phase of the quasi-classical wave function. 

In the relativistic regime, similar to the nonrelativistic case, 
the standard SFA is only exponentially exact since the Coulomb 
field is neglected during ionization, whereas the ITM 129T - 
32 1 can provide also correct preexponential factors. Can the 
quantitatively correct relativistic ionization probabilities be 
derived via the SFA technique accounting Coulomb field effects 
accurately? The relativistic generalized eikonal- Volkov wave 
function (taking also into account quantum recoil) has been 
derived in [66]. The Coulomb corrected SFA based on this 
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wave function has been proposed in [67 1 . However, final results 
have been obtained only in Born approximation, i.e. via an 
expansion of the eikonal wave function with respect to the 
Coulomb field, which, in fact, reduces the transition matrix 
element to the one in the standard second order SFA. 

With this paper we begin a sequel of papers in which 
we develop the relativistic Coulomb-corrected SFA based on 
the Dirac equation, generalizing the nonrelativistic theory of 
l55ll63l and apply it for the calculation of spin-resolved quanti- 
tatively correct ionization probabilities. Rather than the Volkov 
wave function, the eikonal- Volkov wave function is employed 
as final state of the Coulomb corrected SFA. The influence of 
Coulomb potential of the atomic core on the ionized electron 
continuum dynamics is taken into account via the eikonal ap- 
proximation. The latter means that the quasi-classical (WKB) 
approximation is applied for the electron continuum dynamics 
and, additionally, the Coulomb potential is treated perturba- 
tively in the phase of the quasi-classical wave function. The 
formalism is applied for direct ionization of a hydrogen-like 
system in a strong linearly polarized laser field. 

In this first paper of the sequel, we begin with the nonrel- 
ativistic Coulomb-corrected SFA to show in the most simple 
case the scheme of the Coulomb-corrected SFA. The SFA 
formalism is applied to treat the Coulomb field effect of the 
atomic core during ionization systematically and to obtain 
quantitatively correct results which, in particular, for the total 
ionization rate coincide with the PPT result. Two versions of 
the theory based on the velocity and length gauge, respectively, 
are considered. Comparison with the PPT theory is carried out 
and the physical relevance of the two versions is discussed. A 
conclusion is drawn concerning the scheme of the relativistic 
generalization of the Coulomb-corrected SFA. In the second pa- 
per of the sequel, the relativistic Coulomb-corrected SFA will 
be developed, and the next paper in the sequel will be devoted 
to spin effects in relativistic above-threshold ionization. 

The plan of the paper is the following: In section [IT] the 
nonrelativistic Coulomb-corrected SFA in the length gauge 
is considered and differential and total ionization rates for 
hydrogen-like systems are derived. The next section is dedi- 
cated to the Coulomb-corrected SFA in velocity gauge. The 
comparison of the different versions of the Coulomb-corrected 



SFA is carried out in Sec. IV and the conclusion is given in 
Sec. E 



II. NONRELATIVISTIC COULOMB-CORRECTED SFA IN 
THE LENGTH GAUGE 



In this section we show how the nonrelativistic Coulomb- 
corrected SFA in the length gauge is developed. Rather than 
the usual Volkov wave function, it employs the eikonal- Volkov 
wave function to describe the electron continuum dynamics 
accurately, taking into account the Coulomb field effect of the 
atomic core. As we will see in this way the PPT ionization 
rates can be recovered within the SFA formalism. 



A. The standard SFA 

We consider a highly-charged hydrogen-like ion interacting 
with a laser field. The dynamics is governed by the Hamilto- 
nian 

H = H +H int , (1) 

where Hq is the Hamiltonian of the atomic system 

Ho = p 2 /2 + V(r), (2) 

with the atomic potential V(r), the momentum operator p and 
coordinate vector r (atomic units are used throughout). The 
interaction Hamiltonian due to the laser field in length gauge is 



H im (f) = r ■ E(f), 



(3) 



with the laser electric field E(f). The time evolution operator 
U(t, to) of the atom in the laser field can be formulated via the 
Dyson equation 



Jlu 



U(t,t ) = U (t,to) - i dtU(t,t')H in ,(t')Uo(t',to) (4) 



where Uo is the time evolution operator of the atomic system 
without the laser field. The matrix element for a laser induced 
transition from the initial atomic ground state 10(f)) = \0)e' Ipl , 
with the ground state energy —I p , and the ionization potential 
I p = k 2 /2, into a continuum eigenstate of the total system 
l<Ap(0) w i m an asymptotic momentum p is then given by 



M p = -i f 

J -ex 



dt{^ v (t)\H int (t)\4>(t)). 



(5) 



In the SFA, the final continuum state is approximated by a 
Volkov state |i^p (f)), i.e. an eigenstate of a Hamiltonian, where 

In 



(6) 



the electron is only interacting with the laser field ||33 
coordinate space it is given by 



<r| ( A p v (f)) = exp[/5[ ) 0) (r,f)]/(2;r) 



s.3/2 



The function in the exponent S Q °\r,t) - (p + A(f)) • r + 

J t °° dt' (p + A(f')) 2 /2 is the classical action of an electron in 
a laser field in the length gauge. Note that the Volkov wave 
function coincides exactly with the wave function in the zeroth- 
order WKB approximation for the system. The ionization 
matrix element in the SFA yields: 



Mr 



dt(p+A(t)\H int (t)\0)exv[-iS(t)] 



(7) 



with S(t) = jr°° dt'[(p + A(f')) 2 /2 + k 2 /2]. In the adiabatic 
regime, when the laser frequency to is smaller than the ground 
state energy I p and the ponderomotive potential U p - E^/Aa) 2 , 
with the laser field amplitude Eq, the time integration in Eq. 
(|7| can be carried out in good accuracy via the saddle point 
method (SPM), see, e.g., |gjD. This yields 



Mr, 



2n 
i§(t s ) 



-<p + A(r s )|// !nf (f s )|0> exp[-;S (f,)], (8) 
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where t s are the so-called saddle points of the integrable func- 
tion defined by S (t s ) = 0. After a partial integration in Eq. J7J, 
the transition operator in the matrix element can be transformed 
fromH int to V(r) ll69l: 



and the total ionization rate: 



M P = -i V J-^-(v + A(t s 



s )\V(r)\0)exp[-iS(t s )]. (9) 



In the case of a long laser pulse the differential ionization rate 
is expressed via the matrix element as follows 



dw to 



(10) 



where the summation in Eq. ([H} is carried out only over the 
saddle points of one laser period. 



B. SFA for a negative ion 

The calculation of ionization rates is straightforward in the 
case of ionization of a negative ion. The latter can be modeled 
by a zero-range potential V (z) (r) = -(2n / K)5(r)d,-r, with the 
matrix element (p|V <c) |0 (z) ) — - yfi( /(2ti) |70|. In a sinusoidal 
laser field A(f) = (Eo/w) sin(wf) the saddle point equation 
yields: 



sin(wf s ) = - 



Pe 
E /(o 



+ t 




(11) 



with the Keldysh parameter y = koj/Eq, pe = p ■ e, p± = 
|p - (p • e) e|, e = Eo/|E |. In the tunneling regime (y «: 1) the 
saddle points in one laser cycle can be given approximately via 



a perturbative solution of Eq. (Hi with respect to y: 



OJt, 



- arcsin 



Pe 



u>t S 2 — n + arcsin 



Pe 



+ 1 



■A 



Eq/co 



\E(t )\/co 



\E(t )\/oJ ' 



+ 1- 



(12) 



with \E(to)\ = Eq -J 1 - (ojpe/Eo) ■ Inserting the saddle points 
into Eq.( 10 1, yields the differential ionization probability of a 
negative ion 



dw (z) 



d 3 p 2n 2 \E{t )\ 



exp 



3|£(fo) 



(13) 



Since the ratios pe/(Eq/u>) and p x /(Eo/o)) are smaller than 
one in the case of tunnel ionization, we can expand the func- 
tion in the exponent quadratically in terms of momentum and 
neglect the dependence in the preexponential factor. With this 
we arrive at the differential ionization rate: 



dw (z) 



d 3 p 2tt 2 £ ( 



exp 



2^_ 
"3£ 



E P± 3E, 



3 Pe 



(14) 



3£o 



3/2 



n 2k 5 ' 2 



exp 



2^_ 
"3£ ( , 



(15) 



The SFA ionization rates of Eqs. ( |14[ l and ( |15[ ) for a short 
range potential coincide with the ITM result B71II . The physical 
reason is that neglecting the atomic potential after the electron 
is transferred into the continuum, is justified for negative ions. 



C. SFA for a hydrogen-like system 

In the case of atomic ionization, the Coulomb potential of 
the ionic core cannot be neglected in the electron continuum 
dynamics. Therefore, to obtain an accurate ionization rate, 
the wave function of the continuum state \iff p (t)) in the SFA 
ionization amplitude is approximated by the eikonal wave 
function (instead of the usual Volkov function) which accounts 
for the Coulomb field effect of the ionic core Il55ll63l . 

As we noted in the previous section, the Volkov wave func- 
tion is identical to the electron wave function in the laser field 
in the zeroth order WKB-approximation. A systematic im- 
provement of this state compared to the exact continuum state 
can be achieved employing the WKB-approximation for the 
wave function of an electron exposed to the simultaneous ac- 
tion of the laser and the Coulomb field. From the Schrodinger 
equation for an electron in a Coulomb potential V (c) (r) = —x/r 
and a laser field E(f) 



m t ijf= -yAi/r + V (e V + r-E(f)i/', (16) 
the ansatz iff = e' s ^ n yields the following equation 

_ S = <^ + l *> + ,.E + ?f. ,17) 
2 1 2 

Using the WKB-expansion S = So + jS \ + . . ., we obtain the 
equation 

-So=^ + V (c) + r.E( f ), (18) 

So is the classical action of an electron in the laser field and the 
atomic potential. In the eikonal approximation the partial dif- 
ferential equation for So is solved perturbatively in the atomic 
potential V (c) . The zeroth order solution gives the Volkov- 
action 

1 r°° 

S °\r,t) = ($+A(t))-r+~J dt'(p + A(t')) 2 , (19) 

with A(f) = - JT dt'YL(t'), whereas the first order solution 
reads 



J 00 
dt'V (c) (r( f 



))■ 



(20) 



with the trajectory of the electron in the laser field r(f') = r + 
J dt"p(t") and p(f) = p + A(f). The time t can be interpreted 
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as the time and r as the coordinate of the ionization event. Thus, 
the approximate wave function of the electron continuum state 
in the laser and Coulomb field, which is termed as the eikonal- 
Volkov wave function, in the nonrelativistic regime is 



1 



(2tt) 3 / 2 



exp{iSf\r,t) + iS$\r,t)}. 



(21) 



It takes into account the influence of the atomic potential quasi- 
classically up to first order and will be used in the SFA ampli- 
tude of Eq. d3). 

Let us estimate the applicability of the eikonal approxima- 
tion given by the condition Si' <K S . The perturbed action 
can be estimated ~ j V <c, dT ~ J drx/x ~ log(r£ e /r£,) ~ 

log( V-Efl /-Eb) ~ 1, using the potential V (e) ~ k/te (r£ = r • e), 
the initial coordinate before tunneling r^, ~ v c 6t c , the ve- 
locity v c ~ k, the uncertainty of the initial time St c [in 
the latter, we use the time-width of the saddle-point integra- 
tion 5t c ~ 1/ -\[§(ts) ~ 1 / V^£o] and the tunnel exit coor- 
dinate r£ e ~ k 2 /Eq. While the Volkov-action is estimated 
§ ~ p(t c ) 2 t c + I p t c ~ EqtI + I p t c ~ E a /Eo with the tunneling 
time t c ~ yjuj — k/Eq determined by the Keldysh parameter 
and the atomic field E a = k 3 , the eikonal approximation for the 
nonrelativistic ionization problem is valid when 



Eo , 
— <sc 1. 

E n 



(22) 



Note that Eq/E ci < 1/16 in the tunneling ionization regime for 
a hydrogen-like ion. 

To be able to handle the additional term S Q in the 
SFA transition amplitude, we have to make simplifica- 
tions. The time derivative of S ( ^ given by d t S \r,t) m 
_y( c ) + J^°° af?"p(f")), corresponds to the potential energy of 
the ionized electron in the remote future, after it has escaped 
from the bound state. Since the electron left the atomic system 
after ionization and recollision is not considered here, its po- 
tential energy is vanishing for asymptotically large times and 
therefore it is justified to use Si ( r > t) as 0. Consequently, the 
additional term S ( in the exponent of the amplitude has no 
influence on the saddle point equation and leaves the saddle 
points unchanged [73 1, however, it can change the preexponen- 
tial by a factor exp^/S^r, ?,)]. Thus, the Coulomb-corrected 
SFA amplitude of ionization reads: 



M 



/CO 
dt(p + A(t)\H int (t) expf-iS'V f)]|0 (c) > 
CO 

xexp[-/S(f)], (23) 



where |0 (e) ) is the electron bound state in the Coulomb poten- 
tial. The next task is to find an analytic expression for the new 
preexponential factor for times t — t s . Physically, ^^(r, t) 
corresponds to the sum of potential energies the electron pos- 
sesses on its trajectory. When the electron has left the vicinity 
of the atomic core, the potential energy is small and there are 



no further contributions to 5 1 ! (r, t). Since we consider the 
tunneling regime where Eq/oj » k, this situation sets in at a 



very moment of ionization. Thus, it is justified to expand the 
argument in describing the trajectory of the electron, up to 
second order around the saddle point f„ i.e. around the instant 
of ionization: 



r(f') = r + p(f. s )(f' - t s ) - E(f. s )(f' - t s ) 2 /2. 



(24) 



Further, the momentum distribution of the amplitude is dom- 
inated by the exponential function that is located around the 
laser polarization direction, i.e. we can assume in the pre- 
exponential function p = p E e and p(f,) = iVe. Addition- 
ally, it can be argued that the tunnel ionization starts mainly 
in the area around the laser polarization axis r = r^e, i.e. 
at the outskirts of the atom in direction of the laser elec- 
tric field. This typical value for the initial coordinate of 
the trajectory r is justified via the saddle point condition 
for the integral J d 3 r exp[-;p(f s ) ■ r — Kr\ which leads to 
Tj/r? = P(ts)/(ii<). Thus, the integrand in the expression of 



the Coulomb-correction factor of Eq. (20 1 can be simplified: 



1 



l 



\r E + p E (t s )(f - t s ) - E(t s )(f - t s ) 2 /2\ ' 



(25) 



Furthermore, the motion after the electron has left the barrier, 
contributes only as an unimportant phase in the preexponential 



factor in Eq. ( 23 i and the integration limit can be set at the 
tunnel exit: a)to = - arcsin [pe/(Eu/u>)]. With these simplifi- 
cations the integral in Eq. (|20| can be evaluated: 



expps'Vo] 



1 



1 + Vl +4A 
-1 + Vl +4/1 J 

+ O(A), 



(26) 



with the small quantity A - -r ■ E(?j)/2a- which is of the 



order of y/Eo/E a <K 1, see Eq. (22 1. In fact, one can estimate 
A ~ x c Eq/k 2 ~ v a T c Eo/K 2 ~ -\jEo/E a . We underline that in 



all expressions after Eq. (22 1 expansions in this parameter are 
employed. 

We come to the conclusion that in the nonrelativistic regime 
the Coulomb-corrected SFA amplitude differs from the one in 
the standard SFA by the following Coulomb-correction factor: 



Qn 



Ah 



r-E(f s ) 



(27) 



The transition amplitude can then be expressed in a very simple 
form: 



M 



Xco 
dt(p + A(f)|0 (c) > exp {-i§ (f)} ■ (28) 
DO 

This simple form for the ionization amplitude in length-gauge 
Coulomb-corrected SFA is achieved because the Coulomb- 
correction factor Q nr cancels the dipole interaction factor r ■ E 
in the length-gauge matrix element. 

The occurring matrix element is singular at the saddle point: 



<p + A(f)|0^> 



1 



2^k 5 ' 2 



n [k 2 + (p + A(f)) 2 ] 
1 



2 7rE(f s ) 2 (f - t s ) 2 



(29) 
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where in the last step only the leading order term in Eo/E a is 
retained, and the integral in Eq. (j28j) must be calculated via 
the modified SPM (68), taking into account the pole during 
the integration. Compared to the case of a zero-range potential 
this yields a correction factor in the amplitude of 



The corresponding Volkov wave function describing the free 
electron in the laser field in this gauge is 



ifr v (r,t) 



1 



exp[/p • r + iS(t)]. 



(34) 



M ic) 



2^E a 
\E(t )\ 



(30) 



This correction factor is known from ITM |[35ll but appears to 
be reproducible also with the SFA technique. The differential 
ionization rate in the case of a Coulomb potential of the atomic 
core is 



dw {c) 



4 a>K 6 
-r — r exp 

2 p3 r 



"0 



2E a K ^ 
3Eq Eq 



3E; 



-A 



(31) 



and the total ionization rate yields 



4 a/ — exp 



n E, 



1/2 



2^ 
"3£ 



(32) 



These rates are identical to the PPT-ionization rate 
The momentum distribution of the ionized electrons in the 
non-relativistic regime indicates that the emission of electrons 
with a vanishing final momentum is most probable. The lon- 
gitudinal and the transversal widths of the distribution are 
Ay = y/Eo/E a E()/(jL> and Aj_ = y/E{)/E a K, respectively. 

Concluding this section, within the SFA S-matrix formalism 
and employing the eikonal- Volkov wave function for the de- 
scription of the laser-driven electron continuum dynamics dis- 
turbed by the atomic Coulomb potential, as well as neglecting 
recollisions, one can derive quantitatively correct differential as 
well as total ionization rates that coincide with the expressions 
obtained within the PPT quasi-static theory. In the next section 
we apply the Coulomb-corrected SFA formalism in velocity 
gauge. 



III. NONRELATIVISTIC COULOMB-CORRECTED SFA 
IN VELOCITY GAUGE 

It is well known that the SFA is, in general, not gauge- 
invariant and the SFA in different gauges correspond to dif- 
ferent physical approximations. In this section we calculate 
the ionization rate of an hydrogen-like ion using the Coulomb- 
corrected SFA in velocity gauge. Later, we will compare it with 
the results of the PPT theory and the length-gauge Coulomb- 
corrected SFA to answer the question: in which gauge the 
Coulomb-corrected SFA is more relevant for the calculation 
of the ionization rate of an hydrogen-like ion? We will use 
this information in the next paper for the development of the 
relativistic Coulomb-corrected SFA. 

In velocity gauge the Hamiltonian is given by Eq. ([T| with 
the interaction Hamiltonian 



H int (t) = p-A(f)+A(f)72. 



(33) 



In the case of ionization of a negative ion, the ionization ampli- 
tude in the standard SFA in velocity gauge is given by Eq. |9]) 
where the preexponential matrix element is replaced: 



<p + A(f)|V|0> ^<p|V|0>. 



(35) 



Since the matrix element (p|V|0) is constant and does not 
depend on momentum in the case of a short-range potential, it 
is identical to the one in the length gauge. Therefore, the overall 
ionization amplitude for a negative ion is gauge-invariant in 
the standard SFA. 

In the case of a Coulomb-potential as ionic core the situation 
is different. Here the preexponential matrix-element is not 
a constant and the different momentum dependencies could 
lead to a gauge dependence. The Coulomb corrected SFA 
based on the eikonal- Volkov solution can be developed for the 
velocity gauge similar to that in the previous section. The same 
steps lead to the following final expression for the ionization 



amplitude, cf. Eq. (28 i 



dt{p\Q nr [p-A(t) + A(t) 2 /2]\0 (c) ) 

oo 

x exp {-/5(f)}. (36) 

In contrast to the length gauge calculation, the saddle point 
of § lays not on the singularity of the preexponential matrix 
element and the standard saddle point approximation can be 
applied. It yields for the amplitude 



IE a [p E (pE - 2i«) - 2(Pe - !*:)* arctan (^f)] 



V^||£(fo)l 3/2 



x exp 



E 
3/2 



3|£(fo)l 



(37) 



The ionization differential rate in the velocity gauge Coulomb- 
correct SFA reads 



dw (c) 4k 6 oj 



d 3 p 



n 2 E\ 



exp 



2^ 
"3£ 



Eo P± 3El 



p\ 



x 1 + 



4k 2 4k 



Pi 



arctan 



Pe 



2 \ 



1 + 



2k 



4k 2 (p E \ 2 
+ — — arctan — 

Pe U ' 



1 + 



(38) 



The ionization differential rate in the velocity gauge differs 
from that in the length gauge, see Eq. ( 3 1 1, by the expression 



in the curly brackets in Eq. (38 1. To obtain the total ionization 
rate, the p x -integration can be carried out analytically, but 
/^-integration has to be accomplished numerically. 
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FIG. 1: (a) The ratio of the total ionization rate derived in the 
Coulomb-corrected SFA with respect to the PPT rate vs the param- 
eter 6 = y/Eo/E a /y: SFA in length gauge (black, solid), and SFA 
in velocity gauge (red, dashed); (b) The parameter 6 for different 
nuclear charges Z at a fixed angular frequency to = 0.05 a.u. and 
E /E a = 1/25. The laser intensity is / = 5.6 x 10 19 x (Z/10) 6 W/cm 2 . 



In Fig. [T] (a) we compare the total ionization rate calcu- 
lated within the Coulomb-corrected SFA in the length or ve- 
locity gauge with the PPT rate for different values of the pa- 
rameter 6 - ^Eo/Ealy — (Eo/Ea^^ilp/cj). This parame- 
ter arises since the deviation in the two gauges depends on 
the curly bracket that is a function of Pe/k with the typi- 
cal value for the momentum in laser polarization direction 
Pe ~ Ay = y/Eu/E a E()/a>. While the length-gauge result coin- 
cides with the PPT one, the velocity gauge results tends to the 
PPT-rate only in the limits (5 — > l55l and 6 — > oo, deviating 
from the latter at intermediate values of 6. This is evident from 
Eq. {38}, since the curly bracket goes to one in both limits. For 
intermediate values of the parameter 5, the deviation can be 
larger than a factor of 2. Note that in the tunneling regime the 
parameter 5 can vary in the total range of (0, oo). In Fig.[T](b) 
we show the value of 6 for different nuclear charges Z and 
a suboptical angular frequency. It can be seen that for this 
parameter set the value of 5 lays in an area where the results in 
the two gauges differ significantly. 



IV. COMPARISON OF DIFFERENT APPROXIMATIONS 

In the previous sections we have calculated the ionization 
of a hydrogen-like system in a strong linearly polarized laser 
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FIG. 2: The ratio of the total ionization rate derived in the SFA with 
respect to the PPT rate vs the parameter E /E a : (black, solid) in the 
Coulomb-corrected SFA in length gauge, (green, dash-dotted) in the 
Coulomb-corrected SFA in velocity gauge for y = 0.1, (blue, dashed) 
standard SFA with a Coulomb-potential, (red, dotted) standard SFA 
with a zero-range potential. 



field using the Coulomb-corrected SFA in length and velocity 
gauge. In Fig. [2] we compare the total ionization rates in these 
approximations with the PPT ionization rate for different val- 
ues of Eo/E a . For comparison also the rates in the standard 
SFA are presented using a short-range potential and a Coulomb 
potential. All approximations show the same qualitative be- 
havior, but the absolute values of the rates differ significantly. 
The Coulomb-corrected SFA increases the ionization rate by 
several orders of magnitude. This is in accordance with the in- 
tuitive picture that the Coulomb-potential lowers the tunneling 
barrier and therefore facilitates tunneling. Further, it should be 
mentioned that the Coulomb-correction is only depending on 
Eo/E a , but not, e.g., on I p or a>. 

Thus, from the results of this and the previous sections one 
can conclude that the Coulomb-corrected SFA shows a good 
agreement with the PPT theory only in length gauge. This is a 
message that should be taken into account in the generalization 
of the Coulomb-corrected SFA into the relativistic domain. 



V. CONCLUSION 

We have applied the Coulomb-corrected SFA for ionization 
of hydrogen-like systems in a strong linearly polarized laser 
field. The nonrelativistic regime is considered to show how this 
approximation works and how to use the developed procedure 
for a further generalization of the approximation into the rela- 
tivistic domain. The applied Coulomb-corrected strong-field 
approximation incorporates the eikonal-Volkov wave function 
for the description of the electron continuum dynamics. The 
latter is derived in the WKB approximation taking into ac- 
count the Coulomb field of the atomic core perturbatively in 
the phase of the WKB wave function, i.e., in physical terms, 
the disturbance of the electron energy by the Coulomb field is 
assumed to be smaller with respect to the electron energy in 
the laser field. We have derived an analytical expression for 
the ionization amplitude within the Coulomb-corrected SFA 
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in length and velocity gauge. A simple expression for the 
amplitude is obtained when using the length gauge which is 
due to the fact that the Coulomb correction factor (ratio of 
the Coulomb corrected amplitude to the standard SFA one) in 
this gauge cancels the factor of the electric-dipole interaction 
Hamiltonian in the matrix element. Moreover, a Coulomb 
correction factor coinciding with that derived within the PPT 
theory is obtained. The differential and total ionization rates 
are calculated analytically. The calculated total ionization rate 
in length gauge is identical to the PPT-rate, while in the ve- 
locity gauge it can deviate from the PPT result up to a factor 
of 2. Taking into account that the PPT-rate provides a good 
approximation for experimental results, we can conclude that 
the Coulomb-corrected SFA works successfully in the length 



gauge. The SFA in different gauges, in fact, corresponds to 
different partitions of the total Hamiltonian used to develop 
the SFA |72|. Therefore, one can conclude that the relativistic 
generalization of the Coulomb-corrected SFA, which will be 
carried out in the next paper of this sequel, should be based on 
the partition of the total Hamiltonian that in the nonrelativistic 
limit corresponds to the partition of the length gauge SFA. 
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